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uses  wavelength  division  multiplexing  (WDM)  to  allow  multiple  users  on  a  single  fiber 
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expensive.  A  less  costly  alternative  could  be  achieved  if  Hadamard-Walsh  Code  Shift 
Keyed  (CSK)  encoded  signals  were  used.  An  optical  receiver  was  proposed,  using  the 
Mach-Zender  coupler,  capable  of  demodulating  a  Hadamard-Walsh  CSK  encoded  signal. 
Building  on  previous  work  studying  the  characteristics  of  the  Mach-Zender  coupler,  a 
relationship  between  the  probability  of  bit  error  and  signal-to-noise  ratio  (SNR)  was 
developed  for  a  single  user  and  a  DPSK  optical  signal.  This  relationship  was  then  used  to 
develop  an  understanding  of  the  bit-error  rate  to  SNR  relationship  for  a  multiple-user 
CSK  optical  signal.  Using  the  theoretical  performance  as  a  guide,  a  MATLAB  model  was 
then  constructed  to  investigate  the  sensitivity  of  the  receiver  to  non-ideal  components. 
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I.  INTRODUCTION 

The  majority  of  data  transported  through  modern  U.S.  Naval  platforms  is  carried 
on  copper  wire  or  coaxial  cable.  As  the  amount  of  data  being  transferred  dramatically  in- 
creases, the  need  for  a  more  efficient  methods  of  transport  arises.  Given  the  tremendous 
increase  in  bandwidth  of  fiber  optics,  its  use  to  transport  data  in  naval  platforms  is  be- 
coming standard.  Not  only  does  fiber  optic  communication  offer  a  tremendous  band- 
width, it  offers  many  other  desirable  characteristics  that  include  resistance  to  electro- 
magnetic interference  (EMI)  and  an  excellent  signal-to-noise  ratio  that  make  its  use  on 
naval  platforms  extremely  attractive. 

To  design  the  information  flow  in  naval  vessels  of  the  future,  it  is  important  to 
study  methods  of  optical  data  transmission  and  explore  ways  to  improve  upon  and  exploit 
the  tremendous  bandwidth  available.  In  this  thesis  an  alternative  to  Wave  Division  Multi- 
plexing (WDM)  will  be  developed  and  analyzed.  This  alternative  relies  upon  the  use  of  a 
specialized  coupler,  the  Mach-Zender  coupler,  to  demodulate  the  multiplexed  optical  sig- 
nal into  individual  signal  components. 

The  Mach-Zender  coupler  is  composed  of  two  2x2  optical  couplers,  a  delay  line 
and  a  phase  shift  element.  The  signal  enters  the  first  coupler  input  port  (port  1)  as  shown 
in  Figure  1.1  and  is  equally  divided  between  the  output  ports  A  and  B.  Signal  al  is  then 
delayed  one  bit  period  (T)  and  becomes  signal  a[ .  Both  signals  then  enter  the  second 
coupler  are  combined  and  then  equally  divided  to  the  output  ports  A  and  B.  The  signal 
b[  is  then  phase  shifted  -n  12.  radians  giving  the  outputs  of  the  Mach-Zender  coupler  of 
a2  and  b2 . 

Utilizing  the  outputs  of  the  Mach-Zender  coupler,  a  receiver  can  be  constructed  to 
demodulate  a  single  optical  Differential  Phase  Shift  Keyed  (DPSK)  signal  as  shown  in 
Figure  1.2.  The  outputs  of  the  Mach-Zender  coupler  are  followed  by  a  photodetector 
(PD)  to  convert  the  optical  signal  to  an  electrical  current  by  an  amount  proportional  to  the 
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Figure  1 . 1  The  Mach-Zender  Coupler 

responsivity  (9? )  of  the  photodiode.  This  current  signal  is  further  amplified  by  a  transim- 
pedence  amplifier  of  gain  Zq  and  input  to  an  integrator  over  one  period  (7).  The  energy 
of  the  two  channels  is  then  compared  and  a  symbol  ("1"  or  "0")  is  chosen  based  on  the 
channel  with  the  greatest  energy. 
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Figure  1 .2  Demodulation  of  DPSK  Signal  Using  Mach-Zender  Coupler 


By  combining  multiple  couplers  with  varying  delays  (e.g.  T  and  27),  a  Code  Shift  Keyed 
(CSK)  signal  employing  Hadamard- Walsh  coding  could  be  demodulated.  This  orthogo- 
nally coded  signal  would  allow  multiple  stations  to  use  a  single  fiber  optic  line  simulta- 
neously as  is  currently  being  done  using  Wavelength  Division  Multiplexing  (WDM). 
With  WDM,  however,  there  is  a  significant  cost  in  creating  very  sensitive  filters  neces- 
sary to  demodulate  the  signal  [Ref.  1].  Figure  1.3  shows  a  system  using  multiple  Mach- 
Zender  couplers  that  is  capable  of  decoding  four  Hadamard- Walsh  coded  signals. 
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By  increasing  the  number  of  Mach-Zender  coupler  stages,  this  receiver  can  be  con- 
structed to  demultiplex  increasingly  complex  Hadamard- Walsh  coded  signals. 

In  Chapter  II,  the  theory  behind  the  Mach-Zender  coupler  will  be  discussed,  and 
from  this  theory,  a  mathematical  model  of  the  coupler  will  be  developed.  This  model  will 
include  both  time  and  frequency  domains.  Chapter  III  will  use  the  Mach-Zender  mathe- 
matical model  of  Chapter  II  to  develop  a  single  user  DPSK  encoded  optical  receiver. 
From  this  DPSK  receiver  a  probability  of  bit-error  will  be  determined  based  on  the  sig- 
nal-to-noise ratio  at  the  receiver.  Chapter  IV  will  then  propose  a  multi-user  Code  Shift 
Keyed  (CSK)  optical  receiver  and  develop  the  bit  error  rate  for  the  multi-user  receiver. 
Chapter  V  will  then  compare  the  theoretical  results  of  the  bit  error  rate  versus  signal-to- 
noise  ratio  in  the  single  user  receiver  with  those  obtained  by  a  computer  model  of  the  re- 
ceiver. Chapter  V  will  continue  with  the  computer  model  to  conduct  a  sensitivity  analysis 
of  the  single  user  receiver.  Finally,  Chapter  VI  will  discuss  the  conclusion  of  this  thesis 
and  present  some  areas  for  future  research  with  these  optical  receivers. 


II.  DEVELOPMENT  OF  A  MATHEMATICAL  MODEL  FOR  THE  MACH- 

ZENDER  COUPLER 


A. 


OVERVIEW  OF  MACH-ZENDER  COUPLER 


In  the  following  analysis,  a  mathematical  description  of  the  Mach-Zender  coupler, 
as  has  been  previously  developed  at  the  Naval  Postgraduate  School  [Ref.  2],  will  be  re- 
viewed. This  model  will  be  used  to  analyze  the  demodulation  of  an  optical  DPSK  signal 
using  a  single  Mach-Zender  coupler. 

As  previously  described,  the  Mach-Zender  Coupler  consists  of  two  2x2  optical 
couplers,  a  delay  line  of  some  delay  nT  where  n  is  an  integer  and  T  is  the  period  of  the 
signal,  and  a  device  to  provide  a  phase  shift  of  -  n  12  radians.  Port  1  of  the  first  coupler 
is  the  input  and  Port  A  of  the  second  coupler  and  the  output  port  of  the  final  phase  shift 
element  are  the  outputs. 

The  delay  of  the  Mach-Zender  coupler  can  be  accomplished  by  inserting  a  length 
of  fiber  equivalent  to  the  desired  delay.  Heinbaugh  [Ref.  2]  proposed  using  a  2x2  coupler 
with  a  single  input  and  a  single  output  to  obtain  a  7i/2  phase  shift.  This  mathematical 
model  will  use  a  -nl2  phase  shift. 

The  input  signal  is  defined  as  s(t)  in  the  time  domain  and  S(f)  is  its  Fourier  trans- 
form in  the  frequency  domain.  As  the  review  in  this  chapter  will  show,  the  output  in  the 
time  domain  is  s(t-T)-s(t)  for  port  A  of  the  second  coupler  and  s(t-T)+s(t)  for  the  output 
of  the  phase  shift  element.  This  is  as  illustrated  in  Figure  2.1 . 


'A((s(t-T)-s(t)) 


>A((s(t-T)+s(t)) 


Figure  2.1  Signal  Throughput  of  the  Mach-Zender  Coupler 


B.  DEVELOPING  A  MATHEMATICAL  MODEL 

In  this  section  each  element  of  the  Mach-Zender  coupler  will  be  analyzed  mathe- 
matically and  an  expression  relating  to  that  particular  element's  effect  on  s(t)  will  be  de- 
veloped. 

1.  Coupler  Analysis 

The  optical  coupler,  Figure  2.2,  has  four  ports  labeled  1,  2,  A  and  B.  Each  port 
can  be  used  as  an  input  or  an  output . 
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Figure  2.2  Optical  2x2  Coupler 

Let  the  frequency  domain  inputs  of  the  2x2  coupler  of  Figure  2.2  be  described  by 
the  following  vectors  a  and  b  : 


a 


"«l" 

M 

a2 
a3 

b  = 

b2 

_a4_ 
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The  input  vector  a  and  the  output  vector  b  are  related  by  the  scattering  matrix  S 
defined  as  follows  [Ref.  3]: 
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Equation  2.2  can  be  simplified  if  the  matrices  of  Equation  2.4  are  partitioned  as 
shown  below: 
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After  partitioning,  Equation  2.5  can  be  rewritten  as  follows: 
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Observe  that  in  the  ideal  2x2  optical  coupler,  no  part  of  the  input  signal  is  re- 
flected back  to  the  input  ports.  Specifically,  when  the  only  inputs  are  a,  and  a2 ,  the  out- 
puts are  confined  to  by  and  bA  as  in  Equation  2.6b.  That  is  63  and  b4  consist  of  only  in- 
phase  and  /r  /  2  phase-shifted  components  of  ax  and  a2 . 
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Figure  2.3  Single  Input  2x2  Coupler 

To  apply  this  result  to  a  practical  use  of  the  Mach-Zender  coupler,  consider  Figure 
2.3  in  which  one  input  signal  is  applied  to  one  input  port,  and  the  remaining  port  is  termi- 
nated.   Equation  2.6b  now  reduces  to  become: 
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We  see  that  the  input  a,  is  decreased  by  a  factor  of  i/v2  at  ports  A  and  B  and  the  output 
at  port  B  is  shifted  in  phase  by  n  1 2  radians  from  the  input. 


2.    Delay  Element  Analysis 

As  shown  in  Figure  1.1  the  delay  element  represents  a  time  shift  of  T  in  the  time 
domain  which  corresponds  to  multiplication  by  eilTfT  in  the  frequency  domain.  The  input 
and  output  of  the  delay  element  can  be  modeled  in  the  time  domain  as  in  Figure  2.4  with 
the  output  a[  =  axe'lvfT  in  the  frequency  domain. 


«,(')-> 


a,  (/-r)— > 


Figure  2.4  Time  delay  element 


3.    Phase  Shift  Element  Analysis 

The  phase  shift  of  the  Mach-Zender  Coupler  can  be  accomplished  using  a  2x2 
coupler  with  only  one  input  and  one  output.  The  result  would  be  a  signal  shifted  by  n  1 2 
radians.  For  this  model,  a  phase  shift  of  -nil  radians  was  assumed.  The  phase  shift 
can  be  modeled  in  the  frequency  domain  as  in  Figure  2.5. 
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Figure  2.5  Phase  Shift  Element 


C. 


INPUT  SIGNAL  ANALYSIS 


Using  the  preceding  development,  the  function  of  the  Mach-Zender  coupler  can 
now  be  completely  described.  Consider  an  input  signal,  s(t),  applied  to  an  input  port  of 
the  Mach-Zender  coupler.  The  unused  input  port  (2)  is  terminated.  Let  S(f)  be  the  Fourier 
transform  of  s(t).  The  signal  can  be  applied  to  a  Mach-Zender  coupler  as  shown  in  Figure 
2.6. 
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Figure  2.6  Input  Signal  Through  Mach-Zender  Coupler 


Using  Equation  2.7  the  output  of  the  first  coupler  becomes 


jS(f)/y[2 


(2.8) 


After  the  delay  line  the  inputs  to  the  second  coupler  are  obtained  by  applying  the 
delay  illustrated  in  Figure  2.4.  The  inputs  to  the  second  coupler  are: 
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The  outputs  of  the  second  coupler,  b"  and  b\ ,  are  obtained  by  applying  Equation 
2.6b  with  b\  and  b4  as  the  inputs.  The  following  result  is  obtained: 


bi" 

b; 


1 

V2 


1     J 

J     1 


S{f)e-Jlnfr  1 42 
jS(f)l4l 


l/2(S(f)e-j2^-S(f)) 
j/2(S(f)e->2*1  -S(f)) 


(2.10) 


The  addition  of  the  -n  12  phase  shift  gives  the  final  result  in  the  frequency  do- 


main  as: 

~b<\ 

"l 

0 

K 

0 

-J 

l/2(S(f)e-j2*r-S(f)) 
j/2(S(f)e->2*r-S(f)) 


\S(f)e-^-S(f)) 
(Sine-'^+Sif)) 


(2.11) 
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By  taking  the  inverse  Fourier  transform  of  Equation  2.1 1  the  time  domain  result  is 
obtained: 


s{t-T)-s(t) 
s(t  -T)  +  s(t) 


(2.12) 


By  multiplying  the  matrices  used  to  represent  the  phase  shift,  delay  and  two  cou- 
plers (Equation  2.13),  we  can  derive  a  general  equation  (Equation  2.14(a),  2.14(b))  where 
the  output  of  the  Mach-Zender  coupler  can  be  readily  determined  from  the  inputs  a,  and 


a7, 


X" 

"i 

0  " 

1 

"1  j] 

K 

0 

-j_ 

V2 

J  l\ 

Phase 
Shift 


Coupler 


-j2nfT        q 

0         1 
delay 


S 


a, 


a. 


1     J 
j     1. 

Is' 
Coupler 


(2.13) 


K 

l 

"2 

1 

~2 

a, 


a, 


a^e-'2^  +\)  +  a2(e-jlnrr  - 1) 


(2.14a) 


(2.14b) 


D.       SUMMARY 

In  this  chapter  the  optical  2x2  coupler  was  shown  to  combine  the  signals  of  its  in- 
put ports,  split  the  signal,  and  send  the  result  to  the  output  port  such  that  both  signals 
leave  the  coupler  simultaneously.  It  was  determined  that  if  only  one  input  port  was  used 
then  no  part  of  the  input  signal  was  reflected  back.  It  was  also  shown  that  the  outputs 
from  the  Mach-Zender  coupler  are  (S(f)e-JlKfT  -S(f))/2  and  {S(f)e-j2nrr  +  S(f))  12 
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in  the  frequency  domain  or,  equivalently,  (s(t  -T)-  s(t))  1 2  and  (s(t  -T)  +  s(t))  1 2  in 
the  time  domain. 

In  the  next  chapter  a  receiver  will  be  proposed  using  the  Mach-Zender  coupler 
that  is  capable  of  processing  a  DPSK  encoded  signal.  Using  the  results  of  Chapter  II,  the 
signal  will  be  traced  through  the  receiver  and  its  strength  relative  to  the  receiver  noise 
will  be  analyzed.  Finally,  a  probability  of  bit-error,  Pjj,  will  be  developed  based  on  these 
results. 
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III.  DPSK  DEMODULATION  ANALYSIS 

A.  INPUT  DPSK  SIGNAL 

In  Chapter  I,  an  optical  receiver  capable  of  decoding  a  DPSK  encoded  signal  was 
proposed  (Figure  3.1).  Chapter  II  then  developed  a  mathematical  model  (Equation  2.12) 
of  the  Mach-Zender  coupler  to  be  used  in  this  optical  receiver.  With  these  results,  the 
next  step  is  to  analyze  the  input  of  a  DPSK  signal  and  assess  the  receiver  performance 
versus  the  input  signal-to-noise  ratio  (SNR). 


Mach-Zender 
Coupler 

Delay  =  T 

a2-> 

Photo- 
detector 

aj-> 

AMP 
(Z0) 

04-* 

T 

\a4dt 

0 

as-* 

Choose 
Largest 

S(t) 

b2^ 

b4-+ 

"0' 

DPSK 

signal 

or 
'T 

Photo- 
detector 

(S-K) 

&J-> 

AMP 
(Zo) 

T 

\b4dt 

0 

b5^ 

Figure  3.1  Optical  DPSK  Receiver 

A  given  DPSK  signal  will  be  in  one  of  two  forms  (over  a  period  of  2T  where  T  is 
the  period  of  the  input  signal)  depending  on  the  information  bit  being  transmitted.  For  a 
bit  "1"  being  sent  where  0  <  t  <  2T ,  the  signal  will  be: 


Sl (0  =  ±ApT(t  -T)cos(27rfct  +  6>0)  +  ApT(t)cos(2xfct  +  6>0) . 


(3.1) 


For  a  bit  "0"  being  sent  where  0  <  t  <  2T ,  the  signal  will  be: 


S0(t)  =  ±ApT(t  -  T)cos{l7rfct  +  0O)  ±  ApT(t)cos(2nfct  +  <90) 


(3.2) 


where  fc  is  a  multiple  of  1/T and  pT(t  -  iT)  with  i  an  integer  is  defined  as: 
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pT(t-iT)  = 


[1     iT<t<(i  +  l)T 

0     otherwise. 


(3.3) 


Due  to  the  symmetry  of  the  signal,  an  analysis  of  an  input  1  or  0  will  produce  the 
same  results.  To  simplify  the  calculations  an  input  of  0  (Equation  3.2)  will  be  considered. 


B.   RECEIVER  ANALYSIS 

1.    Output  of  Mach-Zender  Coupler 

In  Chapter  II,  the  output  of  the  Mach-Zender  Coupler  was  determined  in  Equation 
2.12.  Applying  Equation  2.12  to  the  input  signal,  given  the  transmission  of  a  0  (Equation 
3.2),  the  following  output  is  obtained  where  T  <  t  <  TT : 


a2(t) 

AC). 


o 

±  A cos(27rfct  +  60)  ±  A cos(2nfct  +  0Q ) 


(3.4) 


Simplifying  this  expression  yields: 


a2{t) 
hit) 


0 

±^cos(2^fcr  +  (90) 


(3.5) 


which  are  the  outputs  of  the  Mach-Zender  coupler  given  an  optical  DPSK  input  signal. 


2.    Photodetector  Analysis 

The  next  stage  of  the  receiver  is  an  optical  photodetector  (PD)  which  converts  the 
optical  signal  into  an  electrical  current  signal.  The  PD  provides  an  output  current  equiva- 
lent to  the  input  peak  power  (P=A^/2)  multiplied  by  the  responsivity  ( 9? )  plus  the  shot 
noise,  r|s(t)  [Ref.  4]. 

Taking  our  input  from  the  Mach-Zender  coupler  (Equation  3.5)  and  applying  it  to 
our  PD  gives  the  following  signal  in  the  receiver  over  the  interval  where  T  <  t  <2T: 
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*3(0 
b3(0 


(3.6) 


It  is  important  to  note  at  this  stage  that  the  shot  noise  produced  by  the  two  PD  de- 
vices ( rjs  (0   and  rjSh  (J) )  are  independent  of  one  another. 

3.    Amplifier  Analysis 

The  next  stage  of  the  receiver  consists  of  a  transimpedance  amplifier  with  gain  Zq. 
Along  with  a  gain  to  the  input  signal,  this  amplifier  will  introduce  additive  white  gaussian 
noise  (AWGN),  rjamp(t),  to  the  signal.  In  order  to  simplify  the  understanding  of  noise  in 
the  receiver  a  summation  of  total  noise  will  be  necessary. 

Both  the  amplifier  noise,  rjamp(t),  and  the  shot  noise,  tjs  (?)  and  rjg  (/) ,  are  mod- 
eled as  AWGN  with  zero  mean.  If  the  total  noise  in  each  individual  channel  is  rj(t)  then: 


r](0=rjamM)  +  T]s(t) 


(3.7) 


Therefore,  the  noise  content  of  each  channel  following  the  amplifier  stage  will  be  rja{t) 
and  Tjh  (?)  with  a  power  spectral  density  (PSD)  of  No-  Again,  the  noise  in  each  channel 
will  be  independent  of  the  noise  in  the  other  channel. 

Given  this  simplifying  assumption,  the  output  of  the  amplifier  is  simply  the  input 
(Equation  3.6)  plus  the  amplifier  noise  ( J]     (t) )  multiplied  by  the  amplifier  gain  (Zq): 


«4(0 
Kit) 


Z0(9U2/2  +  %(0) 


(3.8) 


4.  Integrator  Analysis 

The  next  step  in  the  analysis  is  to  take  the  amplifier  output  (Equation  3.8) 
and  run  it  through  the  integrator.  This  will  give  us: 
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which  will  simplify  to: 


as(t) 
b5(t) 


l[Z071a(t)]dt 

o 
T 

l[(Z0MA2/2)  +  (Z0Tlh(t))]dt 


(3.9) 


a5(t) 

b5(t) 


i 

\[ZQTJa(t)]dt 

o 

r 

Z0%4277  2+J[Z0776(O]^ 


(3.10) 


To  continue  with  the  analysis  of  the  integrator  output  the  statistical  nature  of  the 
signal  at  this  point  will  be  examined.  Signal  a5(t)  is  defined  as  follows: 


a5(t)  =  a5  +  Na(t) 


(3.11) 


where    a5  =  0  andiVa(0  =  \[ZQija(t)]dt .  Since  A^a(/)  consists  of  a  constant  multiplied 

0 

by  zero  mean  AWGN,  its  mean  is  also  zero.  Therefore  a5  represents  the  mean  and  Na  (t) 
represents  the  noise. 


Similarly  if  b5(t)  is  defined  as: 


b5(t)  =  b5  +  Nh(t) 


(3-12) 


where    b5  =  Z0WA2T 1 2  and  Nh  (t)  =  \[Z0r}h (t)]dt ,  then  b5  represents  the  mean  and 

o 
Nh  (t)  represents  the  noise. 

Next,  for  a  gaussian  random  variable  such  as  our  zero  mean  AWGN  defined  by 
Na(t)  and  Nfj(t),  the  variance  (  <j2  )  is  defined  by: 


a*=*Nh2=Z02NJ, 


(3.13) 
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Now  that  the  mean  and  variance  of  a5(t)  and  b5(t)  has  been  determined,  an 
analysis  of  the  decision  process  to  choose  the  signal  sent  can  be  undertaken. 

5.    Decision  Variable  Analysis 

With  the  mean  and  variance  of  the  signals  to  be  input  to  the  decision  element  de- 
termined, an  examination  as  the  probability  of  an  error  in  signal  reception  can  begin.  A 
symbol  "0"  was  assumed  to  be  sent  to  our  device  initially,  therefore  an  error  will  have 
occurred  if  the  symbol  "1"  is  chosen  as  the  output.  This  probability  of  choosing  a  "1" 
when  a  "0"  is  sent,  Pr{l|0},  will  be  equivalent  to  the  magnitude  of  the  signal  in  channel 
a  being  greater  than  the  signal  in  channel  b  given  that  a  "0"  was  sent  or  ?r{a^>b^\0}. 

If  a  new  variable,  Y,  is  defined  where: 

Y  =  b5-a5  (3.14) 

then  it  is  apparent  that: 

Pr{a5  >  b5\0\  =  Pr{7  <  0|0} .  (3.15) 

From  the  analysis  of  the  integrator  output  it  was  determined  that  a5  and  b5  are 
both  independent  gaussian  random  variables,  and  therefore  Y,  as  the  difference  of  a5  and 
b5 ,  is  also  a  gaussian  random  variable.  Consequently  the  mean  of  Y  or  Y  can  be  ob- 
tained as  follows: 

Y  =  b5-a5  =  Z0MA2T/2.  (3.16) 

Similarly  the  variance  of  Y  ( aY2 )  can  also  be  obtained  as  follows: 

aY2=cTNii2+cTNh2=Z02N0T  +  Z02N0T  =  2Z02NJ.  (3.17) 

Pr{l|0}  or  PrJF  <  0|0}  can  be  defined  by  the  probability  density  function  of  Y  as 
follows: 
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Pr{y<0|0}=  \fy{y\0)dy 


(3.18) 


where  the  probability  density  function,  /v(j>|0) ,  for  a  gaussian  random  variable  is  defined 
as: 


•/>|0)  =  ^=— exp 


2.7MJ, 


-iky-Yf) 

2a! 


(3.19) 


Substituting  Equation  3.19  into  Equation  3.18  gives: 


u  u  1 

Pr{r<0|0}=  \fy(y\0)dy  =   h=-exp 


2;r<7„ 


((y-r)2) 

2a2, 


dy 


(3.20) 


(y-Y)  dy 

If  a  new  variable  x  is  defined  such  that  x  = then  dx  =  -  — .  Using  this 

new  variable  to  substitute  into  Equation  3.20  gives: 


1  Ylay  [         2  "1  i  °° 

Pr{r<0|0}=-7=^    Jexp-^-(-o-^)  =  1=    J 

V2^ro-v  _„        L     2  J  V2^  F/ 


exp 


x 

T 


rfr.       (3.21) 


F 

Equation  3.21  is  now  in  the  same  format  as  the  Q  function,  Q{z) ,  where  z  =  — 


Therefore: 


^ 


Pr{l|0}  =  Q  —    =  Q 


( ZMA2TI?} 


VVyJ 


ZoWJ 


( 


=  Q 


1      2  \JL 

2^    pN0; 


(3.22) 
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C.  PROBABILITY  OF  BIT  ERROR 


The  probability  of  bit  error,  Ph ,  for  this  receiver  can  be  determined  as: 


Ph  =  Pr(l)  x  Pr{0|l}  +  Pr(0)  x  Pr{l|o} . 


(3.23) 


If  a  "1"  or  a  "0"  is  equiprobable  (i.e.,  Pr(l)  =  Pr(0) )  and,  if  by  symmetry  Pr{0|l} 
=  Pr{l|o} ,  then  substituting  Equation  3.22  into  Equation  3.23: 


Pb=Q 


I  2  T 


(3.24) 


Further,  if  the  peak  power  (P)  is  defined  as  P-A^/2  then  Equation  3.24  will  simplify  to: 


Ph=Q 


MP. 


T 


IN 


0/ 


(3.25) 


Defining  the  voltage  signal-to-noise  ratio  (SNR)  as: 


Y      1 
SNR  =  —  =  -WA2 


T 


2N0 


MP. 


2N„ 


(3.26) 


then  Pk  becomes: 


Ph  =  Q(SNR) 


(3.27) 


Figure  3.2  shows  this  dependence,  which  illustrates  the  rapid  decline  in  bit  error  rate  as 
the  signal-to-noise  ratio  exceeds  4  dB. 

D.  SUMMARY 

In  this  chapter  a  mathematical  model  for  a  single  user  DPSK  encoded  optical  re- 
ceiver was  developed  using  the  analysis  of  the  Mach-Zender  coupler  from  Chapter  II. 
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0  2 

SNR(dB) 


Figure  3.2  Pb  vs.  SNR 

From  this  model  a  relationship  of  bit-error  rate  versus  signal-to-noise  ratio  was  then  de- 
veloped. 

In  Chapter  IV,  a  mathematical  model  of  a  multi-user  CSK  optical  receiver  will  be 
developed  in  a  very  similar  manner  to  the  development  in  this  chapter  for  a  DPSK  re- 
ceiver. Again  as  in  this  chapter,  that  model  will  be  used  to  determine  the  performance 
(i.e.,  Pb  vs.  SNR)  of  the  CSK  multi-user  receiver. 
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IV.  CSK  DEMODULATION  ANALYSIS 

A.  INTRODUCTION 

The  preceding  chapter  developed  a  probability  of  bit  error  for  a  DPSK  demodu- 
lated signal.  A  DPSK  demodulated  signal  implies  that  only  one  user  can  send  in  informa- 
tion in  the  optical  fiber  at  a  time.  In  order  to  allow  multiple  users  to  utilize  the  fiber  si- 
multaneously, the  signal  demodulation  must  be  modified.  As  was  introduced  in  Chapter  I, 
a  Code  Shift  Keyed  (CSK)  signal  employing  Hadamard-Walsh  coding  can  be  used  to  al- 
low multiple  users  on  a  fiber  simultaneously. 

The  Hadamard-Walsh  code  is  a  straightforward  method  of  creating  multiple  or- 
thogonal signals.  The  Hadamard-Walsh  basic  2x2  matrix  is: 


1      1 
1    -1 


(4.1) 


Observe  that  the  correlation  of  the  first  row  vector  with  the  second  row  vector  is  zero. 
The  matrix  of  Equation  4. 1  can  be  built  upon  to  create  additional  orthogonal  vectors  as 
follows: 


"1 

1 " 

1    1 "  " 

1 

1 

1 

-1 

1  -1 

zl 

1 = 

fl    1 1 

1 

-1 

1 

-1 

1  -1 

(4.2) 


It  is  apparent  that  the  individual  elements  of  the  2x2  Hadamard-Walsh  matrix  are  used  to 
multiply  the  entire  2x2  matrix  and  Equation  4.2  can  be  written  as: 


1111 
1-11-1 
11-1-1 
1-1-1      1 


(4.3) 
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Again  observe  that  each  row  vector  is  orthogonal  (i.e.,  has  a  zero  correlation)  to  all  oth- 
ers. 

The  4x4  matrix  of  Equation  4.3  consists  of  four  orthogonal  vectors  and  can  be  ap- 
plied to  a  system  in  which  two  bits  or  two  users'  information  must  be  decided  simultane- 
ously. It  is  this  matrix  that  will  be  used  to  design  an  optical  receiver  and  develop  the 
probability  of  bit  error  for  two  users  on  the  same  fiber  simultaneously. 

The  design  of  the  receiver  will  be  the  same  as  that  proposed  in  Chapter  I  and  is 
shown  in  Figure  4.1. 

B.  INPUT  CSK  HADAMARD-WALSH  ENCODED  SIGNAL 

A  given  two-user  CSK  Hadamard- Walsh  encoded  signal,  S(t),  will  be  comprised 
of  the  sum  of  two  orthogonal  signals,  one  each  from  channel  A  and  channel  B.  Therefore: 

S(t)  =  SA(t)  +  SB(t).  (4.4) 

S(t)  will  exist  over  a  bit   time  of  T=4T'  where  T"  is  the  period  of  the  input  signal.  Fur- 
ther: 


r = -  (4.5) 

4 


where  T  was  defined  in  chapter  three  as  the  period  of  a  DPSK  signal.  The  possible  signals 
from  channels  A  and  B  where  0  <  t  <  47"  are: 

Sw  (0  =  +Apr  (t  -  37")  cos(2nfct  +  0O)  +  Apr  (t  -  27")  cos(2;zfcf  +  0O ) 

+  Apr  (t  -  7")  cos(2#c/  +  <90 )  +  Apr  (t)  cos(27rfct  +  00 )  (4.6) 

Su  (t)  =  +Apr  (t  -  37")  cos(2nfct  +  <90 )  -  Apr  (t  -  27')  cos{2?rfct  +  0O ) 

+  Apr  (t  -  7")  cos(27ifct  +  6>0 )  -  Apr  (t)  cos(2/z/"cf  +  0O )  (4.7) 

S0B  (0  =  +Apr  (t  -  37")  cos(2?tfct  +  0O)  +  Apr  (t  -  27")  cos(2nfct  +  0O ) 

-  Apr(t  -  T')cos(27rfct  +  0O)-  Apr(t)cos(2tfct  +  0Q)  (4.8) 
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S0A  (0  =  +APr  (t  -  37")  cos(2#f  r  +  0O )  -  Apr  (t  -  27")  cos{2nfct  +  d0 ) 

-  Apr  (t  -  7")  cos(27tfct  +  d0)  +  Apr  (0  cos(27rfct  +  0O )  (4.9) 

where  the  S  subscript  refers  to  the  bit  transmitted  and  the  channel  it  was  transmitted  from, 
respectively.  For  example,  Sja  (t)  refers  to  the  channel  A  component  of  S(t)  in  which  a  1 
was  transmitted  for  channel  A,  and  Sob  (t)  refers  to  the  channel  B  component  of  S(t)  for 
which  a  zero  was  transmitted  for  channel  B.  The  pr{t-iT)  term  is  defined  as  it  was  in 

Chapter  HI  (Equation  3.3).  The  frequency  fc  is  a  multiple  of  7/7'. 

Again  due  to  the  symmetry  of  the  signal  an  analysis  of  all  possible  inputs  is  not 
necessary,  as  each  will  produce  the  same  result.  Also,  as  each  signal  in  channel  A  is  or- 
thogonal to  all  signals  in  channel  B,  there  will  be  no  interference  between  channel  A  and 
channel  B  signals  in  S(t).  To  aid  in  the  development  of  a  probability  of  bit  error,  Pb,  an 
input  of  a  1  in  channel  A  and  a  0  in  channel  B  will  be  considered.  This  input  corresponds 
to  Equations  4.6  and  4.7  respectively  and  the  input  signal  will  be: 

S(t)=S1A(t)+S0B(t).  (4.10) 


C.  RECEIVER  ANALYSIS 


1.    Output  of  First  Stage  Mach-Zender  Coupler 


In  this  multistage  receiver  the  first  stage  is  identical  to  that  analyzed  in  Chapter  III 
(Figure  3.1),  substituting  T  for  T.  Again  as  in  Chapter  HI  the  output  of  the  first  stage  is 
determined  by  applying  Equation  2.12  to  the  input  signal,  Equation  4.9.  Given  a  1  in 
channel  A  and  a  0  in  channel  B  the  following  first  stage  output  is  obtained  where 
7"</<47": 


sA(t) 

SB(t) 


-cos(27tfct  +  d0) 


2Apr  (t  -  37")  +  2Apr  (t  -  2T')  -  2Apr  (t  -  T'j 
+  2  Apr  (t  -  37")  -  2  Apr  (t  -  T) 


(4.11) 


A  detailed  derivation  of  this  output  is  contained  in  Appendix  A. 
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2.     Output  of  Second  Stage  Mach-Zender  Couplers 

The  second-stage  couplers  are  designed  with  a  2T'  delay  variable  and  conse- 
quently the  analysis  of  their  outputs  is  different  than  that  of  the  first  stage.  However,  with 
only  a  slight  modification,  a  variation  of  Equation  2.12  can  be  used  to  describe  their  out- 
puts. Specifically  for  a  delay  of  2T  Equation  2.12  can  be  rewritten  as  follows: 


a{t) 
b{t) 


s(t-2T')-s(t) 
s(t  -  2T')  +  s(t) 


(4.12) 


Applying  Equation  4.12  to  the  input  of  the  second  stage,  Equation  4.11,  the  fol- 
lowing output  is  obtained  where  3T'  <t  <  AT' : 


S'oaM 

0 

s;A(t) 

i 

l/2(4Acos(2jrfct  +  d0)) 

S'0B(t) 

~2 

l/2(-4Acos(27tfct  +  60)) 

s;B(t)_ 

0 

(4.13) 


which  can  be  simplified  to: 


S'oaM 

0 

s;A(t) 

Acos(27Cfct  +  e0)) 

S'0B(?) 

-Acos(27tfct  +  00)) 

S'1B(t)_ 

0 

(4.14) 


These  equations  represents  the  output  at  the  second  stage  Mach-Zender  couplers  given  an 
optical  multiplexed  CSK  Hadamard-Walsh  encoded  signal. 

2.   Photodetector  Analysis 

The  remainder  of  the  receiver  for  each  channel  is  identical  to  the  DPSK  receiver 
analyzed  in  Chapter  HI;  therefore,  the  analysis  will  be  similar.  In  the  next  stage  of  the  re- 
ceiver the  optical  photodetector  (PD)  converts  the  optical  signal  into  an  electrical  signal. 
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■J 

The  PD  provides  an  output  current  equivalent  to  the  input  power  (P=A  /2)  multiplied  by 
the  responsivity,  9^ ,  plus  the  shot  noise,  r|s(t). 

Taking  our  input  from  the  Mach-Zender  couplers  (Equation  4.13)  and  applying  it 
to  our  PD  gives  the  following  signals  in  the  receiver  over  the  interval  37'  <  t  <  AT' : 


s'o'A*) 

i*w 

Suit) 

<K(A2/2)  +  ?7 

S'o'B(t) 

9*(A2/2)  +  7] 

s"B(t)_ 

%.  <*> 

(0 
.(0 


(4.15) 


It  is  important  to  again  note  at  this  stage  that  the  shot  noise  produced  by  the  all  the 
PD  devices  are  independent  of  one  another. 

3.    Amplifier  Analysis 

The  next  stage  of  the  receiver,  as  was  the  case  in  Chapter  in,  consists  of  a  tran- 
simpedance  amplifier  with  gain  Zq.  Along  with  a  gain  to  the  input  signal,  this  amplifier 
will  introduce  AWGN,  T]amp(t),  to  the  signal.  As  was  the  case  in  the  analysis  of  the  DPSK 
receiver  the  total  noise  in  each  channel  will  be  summed.  Since  both  the  amplifier  noise, 
Tlamp(t),  and  the  shot  noise,  T\s(t),  are  AWGN  with  zero  mean,  the  total  noise  in  each  indi- 
vidual channel  is  r\{t)  and  : 


»K0  =  *U(0 +  *?,(') 


(4.16) 


The  noises,  f]oA(t),  i]1A(t),  r\0B(t),  and  7]yB(f ),  all  have  a  power  spectral  density  (PSD)  of  No 
and  the  noise  in  each  channel  is  independent  of  the  noise  in  the  other  channels. 

Given  this  simplifying  assumption  the  output  of  the  amplifier  is  simply  the  input 
(Equation  4.14)  plus  the  amplifier  noise,  riamp(t),  multiplied  by  the  amplifier  gain,  Z0  and 
can  be  written  as: 
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Z(0 


s".'(0 


'OB' 


SlsiO 


ZoW) 
Z0(^2/2  +  7l/4(0) 
Z0(XA2/2  +  T]0B(t)) 

ZoT?iB(t) 


(4.17) 


4.  Integrator  Analysis 


The  next  step  in  the  analysis  is  to  take  the  amplifier  output  (Equation  4.17)  and 
run  it  through  the  integrator.  This  will  give  us: 


which  will  simplify  to: 


S'07(t) 

s,7(0 


s'o7(0 
s[7(t) 
s'Q7(t) 


\[Ztf0A{t)\dt 

iT' 

J[Z0(9W2/2  +  i7M(/))]dlr 

r 

7" 

J[Z0(9L42  /2  +  rf0B(t))]A 

47- 
\[Z0TJXB(t)]dt 


4  7" 


37" 
4  7" 


37" 


3r 


47- 


j[Z0TJ0A(t)]dt 

37" 

47" 

Z0mA2T'/2+  \[Z0TjlA(t)]dt 

37' 

47" 

Z09U2r/2+  J[Z0i70/r(0]A 

47- 

J[Z0i7w(0]* 


3r 


(4.18) 


(4.19) 


To  greatly  simplify  the  analysis  at  this  point,  notice  the  similarity  of  the  integrator 
output  defined  in  Equation  4.19  to  that  of  the  DPSK  receiver  integrator  output  defined  in 
Equation  3.10.    Specifically,  the  information  plus  noise  channels  in  both  cases  (S"A"(t) , 
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S'o'b '(Oand  b5(t))  are  identical  as  are  the  noise-only  channels  (S"B"(t),S'0'A"(t),and  a5(t)). 
Therefore,  the  statistical  nature  of  the  CSK  receiver  integrator  output  can  be  summarized 
as: 

S0A=S0A+N0A  (4.20) 

S1A=S1A  +  NiA  (4.21) 

Sob  =  Sob  +  Hob  (4-22) 

Sim  =  S»  +  Af1B  (4.23) 


where  the  mean  of  each  channel  can  be  described  as  S0A  and  SlB  (which  are  both  equal 
to  zero)  and  SlA  =  S0B  =  Z0%42r72.  Further,  the  variances  of  each  channel  have  identi- 
cal characteristics,  and  as  in  Chapter  DI,  each  channel  represents  a  gaussian  random  vari- 
able. Therefore  the  variance  in  each  channel  can  be  defined  as: 

GN2=Z02Nor.  (4.24) 

Now  that  the  mean  and  variance  for  each  channel  has  been  determined,  an  analy- 
sis of  the  decision  process  can  be  undertaken. 


5.    Decision  Variable  Analysis 

With  the  determination  of  the  mean  and  variance  of  the  signals  to  be  input  to  the 
decision  element,  an  examination  of  the  probability  of  an  error  in  signal  reception  can  be- 
gin. Due  to  the  symmetry  of  the  signal  the  probability  of  a  bit  error  in  channel  A,  Pb{A}, 
is  equal  to  the  probability  of  a  bit  error  in  channel  B,  PbfBJ.  Therefore  an  error  in  channel 
A  will  be  considered  and  the  result  will  be  the  same  for  channel  B.  If  a  symbol  1  was 
sent  in  channel  A,  an  error  will  occur  if  the  symbol  0  is  chosen.  This  probability  of 
choosing  a  0  when  a  1  is  sent  in  channel  A,  PrA{ojl},  will  be  equivalent  to  the  magnitude 

of  the  signal  Soa  being  larger  than  the  signal  Sja-  A  new  variable  is  defined  as  follows: 
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^  _  S\A        Sqa  ■ 


(4.25) 


It  is  apparent  that  to  make  an  incorrect  decision  Y  must  be  less  than  zero  and  therefore: 


Pr,  {0|l}=  Pr,  {S0A  >  SlA  |l}=  Pr,  {y  <  Oil}. 


(4.26) 


In  Chapter  EI  a  probability  of  bit-error  was  developed  (Equations  3.16  to  3.22) 
relating  one  information  plus  noise  channel  to  one  noise-only  channel.  Similarly,  in  this 
CSK  receiver,  channel  A  is  almost  identical  to  the  result  derived  in  Chapter  HI,  with  the 
only  difference  being  that  the  period  in  channel  A  is  now  T  vice  T.  Using  the  result  ob- 
tained in  Equation  3.24  and  substituting  T  for  T  gives: 


Pr„{0ll}=e 


2         V2Nn 


v 


(4.27) 


Now  using  the  definition  of  T',  Equation  4.5,  and  rewriting  in  terms  of  T: 


Pr*W=2 


2         \Ax2Nr 


(4.28) 


D.  PROBABILITY  OF  BIT  ERROR 


The  probability  of  bit  error  for  this  receiver  can  be  determined  as: 


p„=pM=pM- 


(4.29) 


Since  the  probability  of  an  error  in  channel  A  equals  the  probability  of  an  error  in  channel 
B,  Pb  can  be  generalized  by  considering  channel  A  only.  In  channel  A: 


^A}=  PrA(l)xPr,{0|l}+Pr,(0)xPr^|0^ 


(4.30) 
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If  a  1  or  a  0  in  channel  A  is  equiprobable  (i.e.,  PrA(l)=PrA(0)=l/2)  and  if  by  symmetry 
PrA{0\l}=PrA{l\0}  then  by  substituting  Equation  4.28  into  Equation  4.30: 


Pb=PM=Q 


l-XA>    -I- 

2         \4x2Nc 


(4.31) 


If  the  power  (P)  is  defined  as  P=A2/2  then  Equation  4.31  will  simplify  to: 


p>=Q 


(4.32) 


and  if  the  signal-to-noise  ratio  (SNR)  is  defined  as  it  was  in  Equation  3.25,  then  Pb  be- 
comes: 


Pb=Q 


(SNR} 


(4.33) 


This  result  implies  a  3  dB  degradation  in  performance  when  compared  to  the  result  for 
DPSK,  Equation  3.26.  Figure  4.2  shows  a  plot  of  the  theoretical  performance  of  the 
DPSK  and  the  CSK  receivers  and  at  bit-error  rates  below  about  3xlO"3  the  DPSK  receiver 
is  degraded  by  3  dB. 

E.  SUMMARY 

This  chapter  built  upon  the  development  of  the  single-user  DPSK  receiver  from 
Chapter  HI  and  has  shown  that  the  performance  of  a  two-user  CSK  optical  receiver  is  3 
dB  less  than  that  of  a  single  user  DPSK  optical  receiver.  The  next  step  in  the  study  in  the 
use  of  the  Mach-Zender  coupler  in  an  optical  receiver  will  be  to  create  a  computer  model 
of  the  simple  DPSK  receiver  and  obtain  performance  data  from  the  model.  If  the  model 
performance  data  closely  resembles  the  theoretical  data  derived  in  Chapter  HI  then  vari- 
ous components  of  the  model  can  be  altered  to  conduct  a  sensitivity  analysis  of  the  optical 
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Figure  4.2  Pb  vs.  SNR  (theoretical  DPSK  and  CSK  receivers) 

DPSK  receiver.  Chapter  V  will  discuss  the  construction  of  the  computer  model  and  show 
the  results  of  the  sensitivity  analysis. 
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V.  SIMULATION  AND  SENSITIVITY  ANALYSIS 
A.  SIMULATION 

1.  Design  Approach 

The  third  chapter  developed  the  theoretical  probability  of  bit  error  for  a  DPSK 
optical  signal  demodulated  using  the  Mach-Zender  coupler.  If  a  simulation  of  the  optical 
and  electrical  receiver  components  can  be  constructed  using  MATLAB  code  and  the  sta- 
tistical outcomes  of  the  simulated  receiver  at  various  signal-to-noise  ratios  accumulated, 
then  a  comparison  of  the  simulated  performance  to  the  theoretical  will  confirm  the  accu- 
racy of  the  model. 


2.  Mach-Zender  Coupler  Simulation  Construction 

Previous  work  at  the  Naval  Postgraduate  School  [Ref.  2]  has  extensively  studied 
the  physical  characteristics  of  the  Mach-Zender  coupler.  Given  this  understanding  of  the 
coupler,  a  modularized  model  of  the  coupler  was  constructed  in  MATLAB  code.  Appen- 
dix B  contains  a  copy  of  the  code  for  each  module.  Each  element  of  the  coupler  was  mod- 
eled as  an  individual  function  with  varying  operational  characteristics. 

The  optical  2x2  couplers  were  designed  with  a  variable  input  scattering  matrix. 
This  allows  for  varying  amounts  of  coupling  of  the  input  signal(s)  into  the  output  signals. 
In  addition,  a  delay  element  was  constructed  capable  of  creating  any  desired  delay  onto  an 
input  signal. 

The  phase  shift  element  was  modeled  as  another  2x2  coupler  as  was  proposed  in 
the  phase  shift  element  analysis  of  Chapter  n.  A  drawback  of  using  another  2x2  coupler 
to  phase  shift  the  signal  is  the  3  dB  attenuation  of  the  signal  by  the  coupler.  In  order  to 
balance  the  final  output  power  of  both  channels  (£/  and  b"  of  Figure  2.6)  an  attenuation 

element  was  inserted  in  the  non-phase  shifted  output.  As  with  the  other  modules  the 
amount  of  attenuation  was  a  controlled  variable. 
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3.    Electrical  Receiver  Element  Simulation 


The  next  stage  of  the  proposed  DPSK  optical  signal  receiver  converts  the  de- 
modulated optical  output  of  the  Mach-Zender  coupler  into  an  electrical  signal  for  proc- 
essing to  obtain  a  digital  output  signal  as  described  in  Chapter  EH.  Each  electrical  compo- 
nent was  again  designed  as  a  module  with  varying  inputs. 

The  photodetector  can  apply  the  responsivity  (9t )  and  add  additive  white  gaussian 
noise  to  its  output.  The  transimpendance  amplifier  uniformly  amplifies  the  signal  and  also 
adds  additive  white  gaussian  noise.  The  integrator  integrates  the  signal  over  one  period. 
The  relative  strengths  of  the  signal  in  channel  A  and  channel  B  are  then  compared  to 
make  a  decision  of  a  "1"  or  a  "0". 

4.  Statistical  Analysis 

The  MATLAB  model  of  the  receiver  was  repeatedly  run  at  various  SNR  to  obtain 
the  probability  of  bit  error  for  specific  SNR.  The  simulated  results  and  the  theoretical  re- 
sults are  shown  in  Figure  5.1.  The  simulated  results  differed  by  about  4/10  of  a  dB  from 
the  theoretical.  This  indicates  that  the  MATLAB  model  was  accurate. 

With  the  accuracy  of  the  model  confirmed,  the  model  can  be  used  to  vary  the  op- 
eration of  various  components  and  investigate  the  sensitivity  of  receiver  performance  to 
less  than  ideal  components.  Specifically,  the  amount  of  delay  and  the  phase  shift  were 
altered  to  non-ideal  values  and  new  receiver  probability  of  bit  errors  were  compiled. 

B.  SENSITIVITY  ANALYSIS 

1.  Delay  Error 

In  building  the  Mach-Zender  coupler  the  delay  element  was  implemented  by  in- 
serting an  extra  length  of  fiber  equal  in  length  to  one  bit  period.  At  100  Mbps  Heinbaugh 
[Ref.  2]  showed  that  tolerances  of  ±  0.5%  were  achievable.  Inserting  a  delay  error  of  of 
0.5%  into  the  MATLAB  model  showed  no  difference  in  Pb  vs.  SNR.  In  order  to 
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Figure  5.1  Pb  vs.  SNR  (ideal  receiver) 


acquire  measurable  data  delay,  errors  of  5%,  10%  and  20%  were  introduced  into  the 
model.  The  results  are  shown  in  Figure  5.2.  With  5%  delay  error  the  required  SNR  for  a 
specific  Pb  increased  one  half  of  a  dB.  A  10%  delay  error  resulted  in  over  a  3  dB  degra- 
dation in  performance.  At  20%  delay  error  the  receiver  remained  a  random  number  gen- 
erator (i.e,  Pb=0.5)  up  through  a  13  dB  SNR. 
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Figure  5.2  Pb  vs.  SNR  (simulated  delay  error) 


1.    Phase  Shift  Error 

The  phase  shift  error  is  defined  by  the  ability  of  the  2x2  coupler  to  change  the 
phase  of  the  input  signal  while  keeping  attenuation  as  close  to  3  dB  as  possible.  In  the 
simulation,  attenuation  was  maintained  at  3  dB  but  the  effects  of  phase  delay  errors  of 
10°  and  30°  were  investigated.  The  results  are  shown  in  Figure  5.3. 
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Figure  5.3  Pb  vs.  SNR  (phase  shift  error) 


For  a  10°  phase  shift  error  receiver  performance  degradation  was  minimal  and  almost 
zero.  For  a  30°  error  the  receiver  degradation  was  about  2  dB. 


C.  SUMMARY 

This  chapter  reviewed  the  construction  of  a  computer  model  in  MATLAB  code  to 
simulate  the  DPSK  optical  receiver.  The  performance  of  the  model  was  compared  to  that 
predicted  by  the  theoretical  results  obtained  in  Chapter  III  and  the  model  performance 
closely  resembled  theoretical  performance.  Having  established  the  accuracy  of  the  model, 


37 


the  model  was  then  used  to  determine  the  sensitivity  of  the  DPSK  optical  receiver  given 
errors  in  delay  and  phase  shift  elements.  Chapter  VI  will  summarize  the  modeled  results 
and  discuss  areas  for  future  reasearch  regarding  the  use  of  the  Mach-Zender  coupler  in  an 
optical  receiver. 
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VI.  CONCLUSION 

A.  DISCUSSION  OF  RESULTS 

This  thesis  proposed  the  use  of  the  Mach-Zender  coupler  in  a  single-user  optical 
receiver  and  mathematically  illustrated  its  theoretical  implementation.  Next,  building  on 
this  implementation,  the  theory  supporting  the  use  of  the  Mach-Zender  coupler  to  de- 
modulate multiplexed  fiber-optic  signals  was  developed  and  also  illustrated  mathemati- 
cally. From  this  analysis,  we  created  a  computer  simulation  of  the  single-user  optical  re- 
ceiver. Comparing  the  performance  of  the  simulation  with  the  expected  theoretical  per- 
formance established  the  accuracy  of  the  model.  With  an  accurate  computer  simulation 
developed,  the  characteristics  of  the  model  were  then  altered  to  investigate  the  sensitivi- 
ties of  the  optical  receiver  to  non-ideal  components. 

Chapter  IV  showed  theoretically  how  the  optical  receiver  could  be  expanded  to 
demodulate  a  two-user  multiplexed  fiber  optic  signal.  A  two-user  optical  receiver  experi- 
enced a  3-dB  degradation  in  performance  as  compared  to  the  single-user  optical  receiver 
proposed  in  Chapter  HL  The  receiver  of  Chapter  IV  could  be  extended  to  increase  the 
number  of  users  multiplexed  on  a  single  optic-fiber  by  extending  the  number  of  Mach- 
Zender  coupler  stages  with  varying  delays. 

Chapter  V  discussed  the  receiver  simulation  using  a  MATLAB  model,  and  its  per- 
formance closely  followed  theory.  This  indicated  that  the  model  was  accurate  and  that 
adjusting  parameters  of  the  model  to  more  closely  resemble  an  actual  Mach-Zender  cou- 
pler should  give  reasonable  results  in  the  sensitivity  analysis. 

The  sensitivity  analysis  conducted  with  the  model  showed  only  one-half  dB  deg- 
radation from  a  5%  delay  error.  As  was  discussed  in  the  Chapter  V,  a  5%  error  was  well 
above  the  0.5%  error  achieved  in  the  construction  of  the  Mach-Zender  coupler  by  Hein- 
baugh  [Ref.  2].  Therefore,  the  effects  of  delay  error  due  to  construction  on  receiver  per- 
formance should  be  minimal. 

The  sensitivity  analysis  then  varied  the  phase  shift  error  up  to  30°  and  a  degrada- 
tion of  less  than  2  dB  was  experienced.  A  phase  shift  error  of  30°  is  large;  larger  phase 
shift  errors  are  unlikely.  Therefore,  this  receiver  design  does  not  appear  to  be  sensitive  to 
probable  errors  in  phase  shift. 

This  sensitivity  analysis,  conducted  using  the  MATLAB  model,  indicates  that  this 
is  a  realizable  receiver  with  a  strong  potential  of  simultaneously  demodulating  multiple 
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signals  in  one  fiber  optic  line.  If  this  method  proves  to  be  feasible,  it  could  provide  a  cost- 
effective  alternative  to  demodulating  multiplexed  fiber-optic  signals  over  the  current 
technique  using  Wavelength  Division  Multiplexing. 

B.  FUTURE  RESEARCH 

Future  research  on  the  feasibility  of  using  the  Mach-Zender  coupler  should  focus 
on: 

1.  Using  the  existing  Mach-Zender  couplers  to  build  a  DPSK  optical-signal  re- 
ceiver in  order  to  confirm  the  sensitivity  analysis  already  conducted  with  the 
MATLAB  model. 

2.  Constructing  a  CSK  optical-signal  receiver  to  investigate  the  reality  of  de- 
modulating information  from  two  or  more  users  simultaneously  on  one  fiber. 
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APPENDIX  A:  ANALYSIS  OF  OUTPUT  OF  TWO  MACH-ZENDER 
COUPLERS  FOR  CSK  ENCODED  SIGNAL 


1.  SIGNAL  INPUTS 

The  following  two  signals  are  input  to  the  receiver  of  Figure  4. 1  as  S(t): 

SlA  (0  =  +Apr  (t  -  37")  cos(27tfct  +  0O )  -  Apr  (t  -  27")  cos(2^Tcr  +  0O ) 

+  Apr(t-T')cos(2jrfct  +  00)-Apr(t)cos(2jtfct  +  d0)  (A.l) 

SQB  (t)  =  +APr  (t  -  37")  cos(2^fcr  +  60 )  +  Apr  (t  -  27")  cos(2;rfcr  +  60 ) 

-  Apr  (t  -  7")  cos(27tfct  +  60)-  Apr  (t)  cos(2Trfct  +  00 )  (A.2) 

due  to  the  orthogonality  of  the  signals  and  superposition  each  signal  can  be 
considered  separately  and  the  end  results  summed. 

2.  CHANNEL  A  INPUT  ANALYSIS 

First  consider  the  signal  for  channel  A.  To  aid  in  analysis,  cos(27tfct  +  d0)  will  be 

factored  out  and  left  out  of  the  equation  until  the  end  result  is  obtained.  The  output  of 
the  first  coupler  with  a  delay  of  7"  for  the  channel  A  input  will  be  described  as  S%  (t) 
and  Sg  (t) .  The  superscript  refers  to  the  contributing  input  channel  and  the  subscript 
refers  to  the  output  channel.  Therefore  from  7"  <t  <  AT'  the  output  is: 

SAA(t)  =  ^[S]A(t-T')-S]A(t)]  = 

-  [+Apr  (t  -  47")  -  Apr  (t  -  37")  +  Apr  (t  -  27")  -  Apr  (t  -  7") 

-  [+Apr  (t  -  37")  -  Apr  (t  -  IT')  +  Apr  (t  -  7")  -  Apr  (t)]] 
=  -x2x  [-Apr  (t  -  37")  +  Apr  {t  -  27")  -  Apr  (t  -  7")]  (A.3) 
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SAR{t)=X-[S1A(t-T')+S1A(t)}  = 

k+Apr  (t  -AT')-  ApT,  (t  -  37")  +  Apr  (t  -  27")  -  Apr  (/  -  7") 

+  [+Apr  (t  -  37")  -  Apr  (t  -  IT)  +  Apr  (t  -  T)  -  Apr  (/)]] 

=  ix2x[0]  =  0  (A.4) 


Equation  A.4  is  zero  and  can  be  disregarded.  Therefore,  applying  the  result  of 
Equation  A.3  to  the  input  of  the  second  stage  coupler  with  a  delay  of  2T  gives  an  output 
from  3T'<t<  AT  of: 


S^(t)  =  ^[SA(t-T')-SA(t)]  = 

-[-Apr  if  -  ST)  +  Apr  {t  -  AT)  -  Apr  (t  -  37") 

-  \-ApT  (t  -  37")  +  APr  (t  -IT')-  Apr  if  -  7")]] 
=  ^x2x[0]  =  0  (A.5) 


SxV(t)  =  -[SA(t-T')+SA(t)) 


•  A  i 

1 


[-Apr  if  -  57")  +  Apr  (t  -  AT')  -  Apr  if  -  3T') 


+  [-Apr  (t  -  3T')  +  Apr  (t  -  2V)  -  Apr  (t  -  7")] 
=  -  x  2  x  [-Apr  (t  -  37")]  (A.6) 


Equations  A.5  and  A.6  represent  the  output  of  the  second  stage  coupler  for  channel  A. 
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3.  CHANNEL  B  INPUT  ANALYSIS 

Next  consider  the  signal  for  channel  B.  To  aid  in  analysis,  cos{27tfct  +  60)  will  be 

factored  out  and  left  out  of  the  equation  until  the  end  result  is  obtained.  Therefore  at  the 
output  of  the  first  coupler  with  a  delay  of  T'  for  the  channel  B  input  from  T' <t  <  AT'  is: 

SBA(t)  =  ^[SOB(t-T')-SoB(t)]  = 

I[+APr(f-4r)+M.,-3r)-APr((-2rVAPro-r) 

-  [+Apr  (t  -  37")  +  Apr  (t  -  2T')  -  Apr  (t  -  T')  -  Apr  (t)]] 

=  l-x2x[-Apr(t-2T')]  (A.7) 

SBB(t)  =  ^[S1A(t-T')  +  S1A(t)]  = 

+  [+Apr  (t  -  37')  +  Apr  (t  -  2T')  -  Apr  (t  -  T')  -  Apr  (t)]] 
=  -  x  2  x  [+ Apr  (t  -  3T')  -  Apr  (t  -  r )]  (A.8) 


where  the  superscript  refers  to  the  contributing  input  channel  and  the  subscript  refers  to 
the  output  channel.  Equation  A.7  has  only  one  term  and  that  term  and  a  2T'  delayed 
version  of  that  term  will  be  out  of  range  of  3T'  <t<  AT'  and  can  be  disregarded. 
Equation  A.8  can  be  applied  to  the  input  of  the  second  stage  coupler  with  a  delay  of  2T' 
gives  an  output  from  3T'  <  t  <  4T'  of: 

s0BB'(t)  =  ^[sB(t-r)-sB(t)]  = 

^[+APr(t-5T')-APr(t-3T') 

-[+Apr(t-3T')-Apr(t-T')]] 
=  -x2x[-Apr(t-3T')]  (A.9) 
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SlBB'(t)  =  ±[SB(t-T')+SB(t)]  = 

1 


-[+Apr(t-5T')-Apr(t-3T') 


+  [+APr(t-3T')-Apr(t-T')]] 


1 


=  -x2x[0]  =  0. 

2 


(A.  10) 


Equations  A. 9  and  A.  10  represent  the  output  of  the  second  stage  coupler  for  channel  B. 


4.  OUTPUTS 


The  output  of  the  first  coupler  in  equation  4.1 1  is  as  follows: 


sB(t) 


SAA(t)  +  SBA(t) 

s£(t)  +  sBB(t) 


(AM) 


and  consequently  the  output  of  the  second  coupler  in  equation  4.14  is  as  follows: 


X(o" 

Wo 

S'uit) 

^'(0 

S'osit) 

s0V(0 

_S[B(t) 

_O0_ 

(A.12) 


Substituting  Equations  A.5,  A.6,  A.9,  A.  10  into  equation  A.12  will  give  the  result  listed 
in  Equation  4.14. 
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APPENDIX  B:  MATLAB  SIMULATION  PROGRAMS 

This  appendix  contains  the  MATLAB  code  for  the  various  modules  of  the  DPSK 
optical  receiver  model. 


1.  OPTICAL  DPSK  RECEIVER 

%  set  input  paramters 

f    =1e6;  %  frequency 

T    =  1/f;  %  signal  period 

Td  =  T;  %  length  of  delay 

SM  =[1,j;j,1];  %  scattering  matrix 

L    =  3;  %  attenuation  (dB) 

R    =  1 ;  %  responsivity 

No  =  1 ;  %  noise  PSD 

Zo  =  1 ;  %  amplifier  gain 

PWR  =  [8000];  %  power  levels 

t  =  linspace(O.T);  %  time  intervals  for  one  signal 

%  load  number  of  runs 
n=2000000; 


Vo  query  users  for  number  of  runs 

yGn  =  input('enter  number  of  simulations  to  run»  '); 


forp=1:length(PWR) 

A  =  sqrt(PWR(p)); 

%  initialize  data  matrices 
xmtbit  =  zeros(1,n+1); 
rcvbit  =  zeros(1,n); 

%  initial  synchronization  signal  to  receiver 

S1=A*cos(2*pi*f.*t); 

phase  =  1 ; 

%  create  1  st  data  signal 
[S2,xmtbit(1)]  =  onechxmt(f  .phase); 
S2  =  A*S2; 

%  update  phase  setting 
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if  xmtbit(1)==  1 

phase  =  -1*phase; 
end 

for  i  =  1:n 
%  create  1  st  data  signal 
[S3,xmtbit(i+1)]  =  onechxmt(f,phase); 
S3  =  A  *  S3; 

%  update  phase 
ifxmtbit(i+1)==  1 

phase  =  -1*phase; 
end 

%  create  next  input  signal 
S=[S3,S2,S1]; 

%  input  S  to  Mach-Zender  coupler 
[xa,xb]  =  mzcoup(S,T,Td,SM,L); 

%  input  MZ  output  to  RCVR 
[rcvbit(i)]  =  rcvr(xa,xb,R,Zo,No,T); 

%  progress  signal  to  next  period 

51  =  S2; 

52  =  S3; 

end 

%  resize  xmbit  to  remove  last  bit 
xmtbit  =  xmtbit  (1:n); 

%  initialize  error  count 
errcnt  =  0; 

%  compute  total  nr  of  errors 
fori  =  1:n 

errcnt  =  errcnt  +  xor(xmtbit(i),rcvbit(i)); 
end 

Pb(p)  =  (errcnt/n); 

SNR(p)  =  10*log10(R*PWR(p)*(sqrt(T/(2*No)))); 

end 
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%  procedure  to  save  data  to  .mat  file 

clear  P  S; 

load  ideal2dat; 

P=[P,Pb] 

S=[S,SNR] 

clear  L  A  R  No  Zo  SM  t  T  Td  i  n  Pb  SNR  errcnt  S1  S2  S3  rcvbit  xmtbit; 

clear  f  PWR  phase  p  xa  xb; 

save  ideal2dat; 

clear 

%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


2.  ONE  CHANNEL  TRANSMITTER 

function  [S.xmtbit]  =  onechxmt(f,phase) 

%  1  CHANNEL  DPSK  TRANSMITTER  one  bit  output 

%     [S.XMTBIT]  =  ONECHXMT(F) 

%     Creates  a  random  transmit  (XMTBIT)  bit  of  one  or  zero. 

%     Based  on  transmit  bit  it  returns  a  signal 

%     (S)  one  period  in  length  that  is  inphase 

%     if  XMTBIT  is  a  zero  and  1 80  degrees  out  of 

%     phase  if  xmtbit  is  one. 


%  decide  which  bit  transmitter  will  simulate 
%     signal  to  send 

T=1/f;  %  compute  period  of  signal 

t=linspace(0,T);  %  generate  100  time  intervals 

x=  rand  -  0.5; 

if  x  <=  0 
xmtbit=0; 
S=phase*cos(2*pi*f.*t); 

%  compute  S  for  XMTBIT  =  0 
else 
xmtbit=1 ; 
S=-1  *phase*cos(2*pi*f  .*t); 

%  compute  S  for  XMTBIT  =  1 


47 


end 

%  Naval  Postgraduate  School,  ECE 
%  Created  by:  LT  Paul  W.  Wehner 


3.  MACH-ZENDER  COUPLER 


function  [a,b]  =  mzcoupsm(S,T,Td,SM1,SM2,SM3,L) 

%  MACH  ZENDER  COUPLER  with  varied  scattering  matrices  output 

Vo     For  a  given  input  signal,  S  (with  a  period  T, 

Vo     three  periods  in  length),  returns  two  signals, 

Vo     A  and  B  (both  one  period  in  length)  where 

Vo     A  =  S(t-T_actual)+S(t)  and  B=S(t-T_actual)-S(t). 

Vo     T_ACTUAL  indicates  the  actual  length 


°/ 

°/ 

°/ 

°/ 

°/ 

%     (i.e. .seconds)  of  the  coupler  delay 

%     element  and  SM#  indicates  the 

%     scattering  matrix  of  the  individual  #'rd 

%     2x2  optical  couplers.  L  is  the  attenuation 

%     in  channel  a  to  equalizer  signal  strength 

%     lost  in  channel  b  due  to  the  last  phase 

%     shift. 


%  determine  #  discrete  quantities  in 
%    a  signal  period  (T)  of  S 
n  =  length  (S)/3; 


%  demux  signal  S  in  1  st  2x2  coupler 
[a1  ,b1]  =  optcoup2(SM1  ,S); 

%  delay  signal  in  channel  a 
[a2,d]  =  delay(a1,T,Td); 


%  adjust  length  of  channel  a  and 
%     channel  b  signals 
a3  =  a2(d-n+1:d); 
b2  =  b1(d-n+1:d); 


%  2nd  2x2  coupler  output 
[a4,b3]  =  optcoup2(SM2,a3,b2); 

%  phase  shft  element  output 
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[q,b4]  =  optcoup2(SM3,b3); 

%  attenuate  channel  a 
[a5]  =  atten(a4,L); 

%  define  MZ  coupler  output 
a  =  real(a5); 
b  =  real(b4); 

%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


4.  OPTICAL  COUPLER 


function  [a,b]  =  optcoup2(SM,S1,S2) 

%  optical  2x2  coupler 

%     OPTCOUP2(SM,S1  ,S2) 

°/ 

°/ 
°/ 

o, 


o, 


Vo  For  a  given  input  signal,  S1 ,  and  a  possible 

Vo  second  signal  S2,  provides  two  outputs, 

Vo  A  and  B.  Both  outputs  are  attenuated  by  a 

Vo  factor  of  3dB  and  multiplied  by  the  scattering 

Vo  matrix,  SM.  SM  must  be  [2x2]. 


%  determime  if  second  signal  S2  was  input 
if  nargin  ==  2 

S2=zeros(1,length(S1)); 
elseif  length(S1)  ~=  length(S2) 

disp('ERROR  -  Input  arrays  (S1  and  S2)  must') 

disp('         be  of  the  same  length1) 

a=0; 

b=0; 

return; 
end 

%  define  input  matrix 
y  =  [S1;S2]; 

%  multiply  by  scattering  matrix 
x  =  SM*y; 

%  attenuate  output  by  3dB 
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x  =  x./sqrt(2); 

%  define  individual  outputs 
a  =  x(1,:); 
b  =  x(2,:); 


%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


5.  DELAY  ELEMENT 


function  [a,d]  =  delay  (a,T,Td) 
%  DELAY  ELEMENT 

%     DELAY(A,T,Td),  for  a  given  signal  A  (A  must 
%     be  three  periods  in  length  with  a  period  of 
%     T  seconds)  outputs  a  new  signal  A  delayed  Td 
%     seconds.  Indicates  at  which  element 
%     value,  N,  in  a  the  delayed  vector  begins 
Vo     where  the  vector  is  listed  right  to  left 


o. 


%  create  working  copy  of  a 
b  =  a; 

%  determine  #  discrete  quantities  of  one 
%    period  of  a 
i  =  length(a)/3; 

%  compute  delay  as  percentage  of  actual  delay 
perc_Td  =  Td/T; 


o, 


Vo  equate  percent  delay  to  #  of  discrete 
Vo     quantities  of  a 
j  =  round(i  *  perc_Td);     %  quantized  delay 
k  =  (3*i)  -  j;  %  new  start  of  a 


%  zero  out  delayed  discrete  elements  of  a 
forl  =  (3*i):-1:(k+1) 

a(l)=0; 
end 
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m  =  3*i; 

%  overwrite  a  with  delayed  version  of  itself 
for  n=  k:-1:(2*i-j+1) 

a(n)=b(m); 

m=m-1; 
end 

d=k; 

%  Naval  Postgraduate  School,  ECE 
%  Created  by:  LT  Paul  W.  Wehner 


6.  ATTENUATOR 


function  [a]  =  atten(a,L) 

%  ATTENUATOR 

%     ATTEN(A,L) 

Vo     For  given  input  signal,  A,  performs 

%     attenuation  on  signal  equal  to 

%     L  dB. 


o. 


a  =  a/(sqrt(10A(L/10))); 

%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


7.    ELECTRICAL  CIRCUIT  ELEMENTS 


function  [rcvbit]  =  rcvr(xa,xb,R,Zo,No,T) 

%  OPTICAL  DPSK  SIGNAL  RECEIVER  OUTPUT 

%     [RCVBIT]=CRCVR(XA,XB,R,ZO,NO,T) 

%     For  given  input  signals,  XA  and  XB 

%     this  applies  photodetection 

%     proportional  to  the  responsivity  (R)  to 

Vo     convert  the  optical  signal  to  electrical, 

Vo     amplifies  the  electrical  signal 

Vo     with  a  gain  of  (Zo)  and  then 


°/ 

°/ 

°/ 

%     integrates  the  signal  over  a 

%     single  time  period  (T)  to  obtain  signal 
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%  energy  in  each  branch.  A  comparison  of 

%  two  branches  is  then  made  to 

%  determine  the  transmitted  bit.  Noise  is 

%  added  in  the  photodetector  and 

%  the  amplifier  with  a  power  spectral 

%  density  of  No/2. 


%  inut  optical  signal  to  photodetector 
xa1  =  pd(xa,R,No,T); 
xb1  =  pd(xb,R,No,T); 

%  input  signal  to  transimpedance  amplifier 
xa2  =  tranamp(xa1  ,Zo,No,T); 
xb2  =  tranamp(xb1  ,Zo,No,T); 

%  integrate  both  channels  over  T 
xa3  =  intg(xa2,T); 
xb3  =  intg(xb2,T); 

%  compare  channels  to  decide  received  bit 
if  xa3  >  xb3 

rcvbit  =  1 ; 
else 

rcvbit  =  0; 
end 

%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


8.    PHOTODETECTOR 

function  [a]  =  pd(x,R,No,T) 

%  PHOTODETECTOR 

%     PD(x,R,No,T) 

%     For  given  input  signal,  X,  performs 

%     photdetection  on  signal  power  proportional 

%     to  the  responsivity,  R,  in  amps/watt. 

%     Assume  signal  is  output  over  a  1  ohm 

%     load  such  that  amps  correlate  directly  to 

%     volts.  Adds  photodetection  shot  noise  to 

%     output  signal  A  with  PSD(No/2)  as  defined 

%     from  user  input. 
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%  Create  shot  noise 
s  =  randn(1  ,length(x)); 

%  adjust  shot  noise  PSD 
amp  =  sqrt(No/T); 
s  =  amp  *  s; 

%  compute  signal  power 
Px  =  (x.A2)./2; 

%  multiply  input  power  by  responsivity 
%  to  get  output  in  amps 
Ax  =  R  *  Px; 

%  add  shot  noise  to  create  output  signal 
a  =  Ax  +  s; 


Vo  Naval  Postgraduate  School 

Vo  Created  by:  LT  Paul  W.  Wehner 


o 


9.    TRANSIMPEDANCE  AMPLIFIER 


function  [a]  =  tranamp(x,Zo,No,T) 

%  TRANSIMPEDANCE  AMPLIFIER 

%     TRANAMP(x,Zo,No,T) 

Vo     For  given  input  signal,  X,  performs 

%     uniform  amplification  on  signal 

%     with  gain  Zo.  Adds  noise  to  signal  prior 

%    to  amplification  with  PSD  No/2. 

Vo     Assume  signal  is  output  over  a  1  ohm 

Vo     load  such  that  amps  correlate  directly  to 

Vo     volts.  Adds  photodetection  shot  noise  to 

Vo     output  signal  A  with  PSD(No)  as  defined 

Vo     from  user  input. 


%  Create  amp  noise 
s  =  randn(1  ,length(x)); 

%  adjust  amp  noise 
amp  =  sqrt(No/T); 
s  =  amp  *  s; 
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%  add  amp  noise  to  input  signal 
a  =  x  +  s; 

%  adjust  output  by  gain,  Zo 
a  =  Zo  *  a; 


%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


10.  INTEGRATOR 


function  [a]  =  intg(x.T) 
%  INTEGRATOR 
%     INTG(x.T) 

%     For  given  input  signal,  X,  performs 
%     trapazoidal  integration  on  signal 
%     over  one  time  period  T.  Assumes 
%     elements  of  X  are  equally  spaced 
%     from  0  to  T. 


a  =  (T/length(x))  *trapz(x); 


%  Naval  Postgraduate  School 

%  Created  by:  LT  Paul  W.  Wehner 


11.  Q  FUNCTION 

function  [Q,P]  =  qfn(x) 

%  This  function  computes  the  Gaussian  Q-function 

%  using  the  rational  approximation  26.2.17  of 

%     Abromowitz  and  Stegun.  For  sufficiently  large 

%     arguments,  the  asymptotic  expansion  26.2. 12  is  used. 

%     The  rational  approximation  is  accurate  to  within 

%     7.5e-8.  Q  is  the  value  of  the  Q-function  (integral  of 

%     Gaussian  pdf)  of  zero  mean  and  unit  variance  from  the 

%     input  quantity  to  infinity,  and  P  =  1  -  Q  is  the 

%     cumlative  distribution  function.  The  call  is 
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%  [Q,P]=qfh(x). 

% 

%  R.E.  Ziemer 

%  4/17/94 

b(l)=0.3 1938 153; 
b(2)=-0.356563782; 
b(3)=l. 78 1477937; 
b(4)=-l. 821255978; 
b(5)=  1.330274429; 
p=0.2316419; 
y=abs(x); 
T=l/(l+p*y); 
TT=zeros(l,5); 
for  r=  1:5 

TT(r)=TAr; 
end 

Z=exp(-yA2/2)/sqrt(2*pi); 
if  y<4 

Ql=Z*b*TT; 
else 

Ql=(Z/y)*(l-l/yA2+3/yA4); 
end 
ifx>0 

Q=Q1; 
else 

Q=l-Ql; 

end 

P=l-Q; 
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